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Neutrino-Induced Production of 9 Be in Core-Collapse Supernovae 
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We present two new mechanisms for ^-induced production of 9 Be in core-collapse supernovae. 
In both cases, v reactions in He shells produce 3 H, which then makes 7 Li through 4 He( 3 H, 7) 7 Li. 
For progenitors of ~8 Mq , the shocked He-shell material rapidly expands, allowing production of 
9 Be through 7 Li( 3 H, no) 9 Be. For ultra-metal-poor progenitors of ~11— 15 Mq, neutrons generated 
by v reactions in He shells lead to both 7 Li(n, 7) 8 Li(n, 7) 9 Li(e~i/ e ) 9 Be and a rapid neutron-capture 
process, depending on details of v spectra, flavor oscillations, and the explosion. We discuss the 
associated production of 7 Li and 11 B, related mechanisms for 9 Be production that might operate in 
other sites, and consistency of the identified mechanisms with the observational data. 
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It was argued four decades ago that interactions be- 
tween Galactic cosmic rays (GCRs) and nuclei in the in- 
terstellar medium (ISM) could approximately account for 
the abundances of 6 Li, 7 Li, 9 Be, 10 B, and n B observed 
in the present Galaxy Among these isotopes, 9 Be has 
been regarded as special. While big bang nucleosynthesis 
(BBN) produced an initial abundance of 7 Li Q and the 
v process in core-collapse supernovae (CCSNe) may ac- 
count for much of the Galaxy's current inventory of 7 Li 
and 11 B and some fraction of its 10 B 3], it is widely ac- 
cepted that 9 Be isproduced almost exclusively in the ISM 
by GCRs (e.g., [4]). In this regard, high- velocity ejecta 
from very energetic CCSNe can be considered similar to 
GCRs, though these events may be too rare to gener- 
ate significant amounts of 9 Be in the ISM [5(. Yet recent 
observations (e.g., [6j) show that there is a linear correla- 
tion between log(Be/H) and log(E/H) with a slope ~0.9- 
1 over ~3dex, where E stands for O, Mg, Ti, and Fe, all 
of which are major products of CCSNe. This correlation 
appears more consistent with primary CCSN production 
of Be than with production by GCRs only (cf. [4{). 

In this letter we present two new mechanisms for 
CCSN production of 9 Be that do not involve interac- 
tions with the ISM. The mechanisms are triggered by in- 
elastic interactions of vs, from the cooling proto-neutron 
star (PNS) with nuclei in the He zone of the CCSN 
progenitor. We calculate this nucleosynthesis with the 
hydrodynamic code KEPLER [7] using the most re- 
cent models of ultra-metal-poor massive progenitors also 
evolved with this code. For a progenitor of 8.1 Mq, the 
He shell, initially at a radius r ~ 10 9 cm, is exposed 
to an intense flux of us during its expansion following 
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the passage of the shock. Production of 9 Be occurs 
through 4 He( 3 H,7) 7 Li( 3 H,ri ) 9 Be with 3 H made by v 
reactions on 4 He. For progenitors of 11 and 15 Mq, it 
has been shown recently Q that the charged-current re- 
action 4 He(f/ e , e + n) 3 H in outer He shells at r ~ 10 10 cm 
produces neutron densities large enough to drive an r- 
process, the mechanism by which heavy nuclei are synthe- 
sized through rapid neutron capture on Fe-group seeds, 
given a favorable effective v e spectrum (see below). We 
show that a "mini-r process" occurs simultaneously to 
yield a correlated production of 9 Be from 7 Li. We explore 
the sensitivities of the above two mechanisms of 9 Be pro- 
duction, as well as the associated i/-process production 
of 7 Li and n B, to v emission spectra, flavor oscillations, 
and the explosion energy. We also discuss similar mech- 
anisms that may produce Be at higher metallicities and 
implications for chemical evolution of Be and associated 
elements. 

In the updated version of KEPLER, a full reaction 
network is used to evolve the nuclear composition of a 
massive star throughout its lifetime and to follow the 
nucleosynthesis that accompanies the explosion. This in- 
cludes the nucleosynthesis associated with shock heat- 
ing of the star's mantle and i/-process nucleosynthesis 
associated with a v burst carrying ~300 B ( "Bethe" ; 
IB = 10 51 ergs). All v reactions on 4 He are included 
as in Ref . Q . The KEPLER progenitors employed here, 
denoted by u8.1, ull, and ul5, have initial metallici- 
ties (total mass fraction of elements heavier than He) 
of Z = 1CT 4 Z and masses of 8.1, 11, and 15 M Q 0]. 
These models are very similar to earlier ones evolved with 
a more limited reaction network 10]. The three selected 
stars develop Fe cores by the end of their evolutions. We 
simulate an explosion by driving a piston into the col- 
lapsing progenitor and following the propagation of the 
resulting shock wave. We assume that the PNS cools by 
v emission according to L v (t) = L u (0) exp(— t/r v ), with 
an initial luminosity per species of L V (Q) = 16.7 B/s and 
time constant r„ = 3 s, where t = marks the launching 
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of the piston. The v spectra are approximated as Fermi- 
Dirac distributions with zero chemical potential and fixed 
temperatures T Ve , T Pcl and T„ x = (x = //, r). 

A large set of calculations was performed to assess the 
sensitivity of the nucleosynthesis to the explosion energy, 
v spectra, and flavor oscillations. Numerical CCSN sim- 
ulations using progenitors similar to u8.1 produce weak 
explosions with energies of -E CX pi ~ 0.1 B [11], while ob- 
servations suggest -E 0X pi ~ 1 B for more massive pro- 
genitors of ~13-20M Q (Fig. 1 of [Uj]). Below we dis- 
cuss results for -E oxp i ~ 0.06-0.3 B (u8.1 models) and 
0.1-1 B (ull/ul5 models). Two sets of v temperatures 
were used: (T Ue , T 9b , T Vm ) = (4, 5.33, 8) MeV (H) and 
(3,4,6) MeV (S), which represent the harder and softer 
spectra obtained from earlier (e.g., [13]) and more re- 
cent [14] v transport calculations, respectively. For an 
inverted v mass hierarchy, v e O v x oscillations occur be- 
fore vs reach the He shell, which greatly increases the 
rate of 4 He(P e , e + n) 3 H and hence, the neutron density in 
the outer He shell Q. We will explore the case of full 
v e O v x interconversion. We label the nucleosynthesis 
calculations by the progenitor model, the v spectra, and 
the explosion energy in units of B, with, e.g., u8.1H.l 
indicating progenitor model u8.1, the harder v spectra 
H, and E cyj? \ — 0.1 B. Calculations including v e <-> v x 
oscillations are denoted by a bar above the H or S. The 
abundances of 28 Si and 32 S in the He shells of models 
ull and ul5 are ~10-30 times larger than those found 
in the older models of Ref. [To[. Consequently, we also 
considered modified models ull* and ul5* in which the 
He-shell abundances of 28 Si and 32 S were reduced to their 
former values. Representative results for the total mass 
yields of 9 Be, 7 Li, n B, and Fe are given in Table |TJ 



TABLE I: Yields of 9 Be, 7 Li, n B, and Fe in M ( 



model 


9 Be 




7 Li 




n B 




Fe 




u8.lH.06 


1.57(- 


-10) 


2.77(- 


-7) 


1.47( 


-7) 


1.89( 


-3) 


U8.1H.1 


1.97(- 


-10) 


2.97(- 


-7) 


1.35( 


-7) 


1.75( 


-3) 


u8.1H.l 


1.20(- 


-10) 


2.47(- 


-7) 


1.34( 


-7) 


1.79( 


-3) 


U8.1S.1 


5.02(- 


-11) 


1.11(- 


-7) 


5.67( 


-8) 


1.79( 


-3) 


U8.1S.1 


2.55(- 


-11) 


8.19(- 


-8) 


5.05( 


-8) 


1.80( 


-3) 


u8.1H.3 


2.56(- 


-10) 


3.03(- 


-7) 


1.06( 


-7) 


1.45( 


-3) 


ullH.l 


1.43(- 


-9) 


2.0-2.3(- 


-7) 


2.2-8.7( 


-7) 


< 7.73( 


-2) 


ull*H.l 


9.14(- 


-9) 


1.5-1.9(- 


-7) 2.6-9.5( 


-7) 


< 7.68( 


-2) 


ull*H.3 


9.81(- 


-10) 


3.26(- 


-7) 


1.09( 


-6) 


< 8.75( 


-2) 


ul5H.l 


< 5.20(- 


-10) 


< 3.33(- 


-7) 


< 1.34( 


-6) 


< 4.42( 


-2) 


ul5*H.l 


< 2.92(- 


-9) 


< 3.15(- 


-7) 


< 1.36( 


-6) 


< 4.42! 


-2) 


ul5*H.3 


7.21 (- 


-10) 


1.69(- 


-7) 


9.50( 


-7) 


< 4.621 


-2) 



Note: X(Y) = X xW 



Progenitors of ^8 M Q have a steeply-falling density 
profile outside the core and He-shell radii r ~10 9 cm, in 
contrast to ^10 10 cm for progenitors of > 11 Mq. We 
use Zone 95 in u8.1H.l to illustrate 9 Be production in 
such progenitors. Prior to shock arrival, the radius, tem- 
perature, and density of this zone are 1.58 x 10 9 cm, 
2.21 x 10 8 K, and 279 g/cm 3 , respectively. The three most 



abundant nuclei are 4 He, 12 C, and 16 with initial mass 
fractions of 0.948, 0.043, and 0.009, respectively. The 
time evolution of the number fraction Yi for various nu- 
clei is shown in Fig. [T] Upon being shocked at t ~ 0.7 s, 
Zone 95 reaches a peak temperature of ^8 x 10 8 K, so 
that any 9 Be produced previously is burned up. By 
t ~ 5 s the shocked material has expanded and cooled to 
^2 x 10 8 K, effectively turning off the principal destruc- 
tion reactions 9 Be(p, 4 He) 6 Li and 9 Be(p, c?)2 4 He. Yet 
as the material is still close to the PNS and the time 
still relatively early in units of t v , the flux of vs can 
efficiently drive the breakup reactions 4 He(^, i/n) 3 Hc, 
4 He(i/, i/p) 3 H, and 4 He(P e , e+n) 3 U. Production of 9 Be 
occurs through 4 He( 3 H, 7 ) 7 Li( 3 H, n ) 9 Be. Note that 9 Be 
must be produced in the ground state (hence no) because 
all of its excited states are unstable to breakup. We took 
the rate for 7 Li( 3 H, no) 9 Be from Ref. [151 ] . 




t(s) 



FIG. 1: Time evolution of the number fraction Yi for various 
nuclei in Zone 95 of u8.1H.l. 

The yield of 9 Be decreases by a factor of ~4 from 
u8.1H.l to uS.lS.l (Table This reflects the high 
thresholds of the breakup reactions 4 He(z^, z/p) 3 H and 
4 He(P e , e + n) 3 H, and consequently, their sensitivity to 
the high-energy tails of the v spectra. The presence of 
neutral-current production of 3 H reduces the sensitivity 
to v e O v x oscillations: other factors being the same, 
oscillations increase the yield by up to a factor of ^2. 
We also explored different explosion energies based on 
current_CCSN models [Il|. The change from u8.1H.l 
to u8.lH.06 (u8.1H3) produces a 20% decrease (30% in- 
crease) in the 9 Be yield primarily through the expansion 
rate of the shocked He-shell material: a higher E c ^ p \ re- 
sults in faster cooling to ~2 x 10 8 K and hence, more 
9 Be. The 9 Be production in progenitors like u8.1 depends 
on the sharply-falling density structure outside the core: 
the He shell is shocked, expands, and cools on timescales 
comparable to r„, and thus before the v flux has dimin- 
ished significantly. In such progenitors the yield is essen- 
tially independent of the initial metallicity. 

We choose ull*H.l to illustrate the other new mecha- 
nism for ^-induced production of 9 Be. This model is es- 
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sentially the same as that used in Ref. Q to demonstrate 
a i/-induced r-process in the He shell at low metallicities, 
a scenario first proposed in Ref. [l6| and later explored 
in Refs. 0, [T3]- The importance of 4 He(£/ e , e+n) 3 H to 
neutron production and the crucial effect of v e <-> v x os- 
cillations were pointed out in Ref. [8] . Prior to shock ar- 
rival, the radius, temperature, and density of Zone 223 in 
ull*H.l are 1.10 x 10 10 cm, 8.49 x 10 7 K, and 50 g/cm 3 , 
respectively. This zone is composed predominantly of 
4 He, with 12 C at an initial mass fraction of ~10~ 5 the 
second most abundant species. The initial mass fractions 
of 28 Si, 32 S, and 56 Fe are ~10~ 8 , 4 x 10~ 9 , and 4 x 10~ 8 , 
respectively. The time evolution of Yi for various nuclei 
is shown in Fig. [2] 

Neutrons produced mainly by He(S' e , e + n) 3 H are cap- 
tured by Fe-group nuclei to drive an r-process in Zone 
223 0]. An analogous "mini-r process" occurs to pro- 
duce 9 Be through 7 Li(ri, 7) 8 Li(n, 7) 9 Li(e _ j> e ) 9 Be. The 
9 Be yield is limited by the short 838 ms half-life of 8 Li, 
which j3 decays through the 3.0 MeV resonance in 8 Be 
to 4 He + 4 He, and by the 49.5% branching ratio for 9 Li 
to decay to particle-unstable excited states in 9 Be. The 
above reaction chain operates for ^20 s prior to shock 
arrival, at which time the temperature and density jump 
to ^2 x 10 8 K and ^220 g/cm 3 , respectively, and a short 
burst of neutrons is released by 8 Li( 4 He, n) 11 B (Fig. [5]). 
Over the few seconds, when the postshock temperature 
stays at ^2 x 10 s K, incomplete destruction of 9 Be oc- 
curs through 9 Be(p, 4 He) 6 Li and 9 Be(p, d)2 4 He, due to 
the absorption of protons from 4 He(j/, v'p) 3 H (Fig. [5]). 
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FIG. 2: Time evolution of the number fraction Yi for various 
nuclei in Zone 223 of ull*H.l. 

In contrast to the u8.1 models, the mechanism for the 
ull/ul5 models produces 9 Be prior to arrival of the shock 
and thus is sensitive to the explosion energy (Table A} . 
The higher postshock temperatures found in more en- 
ergetic explosions greatly enhance 9 Be destruction. For 
Ecxpi ^ 1 B, essentially all of the 9 Bc made in the He shell 
prior to shock arrival is destroyed: only that produced in 
the oxygen shell through 12 C(v, z/) 3 He + 9 Be [3] survives. 
The absence of protons and 4 He in the oxygen shell elim- 
inates light-particle reactions through which 9 Be can be 



destroyed. 

The 9 Be yields for the ull/ul5 models are also sen- 
sitive to the composition of the He shell (Table IJ). 
The abundances of the hydrostatic-burning products 28 Si 
and 32 S are ~10-30 times higher in ull/ul5 than in 
ull*/ul5*. As 28 Si and 32 S are neutron sinks, Y n is con- 
sequently reduced in ull/ul5. The effects on 9 Be yields 
are quadratic in Y n , as two neutrons are captured to pro- 
duce 9 Be. Increasing the initial metallicity also reduces 
Y n due to neutron sinks such as 56 Fe. He-shell produc- 
tion of 9 Be dominates oxygen-shell production only for 
initial metallicities of Z < 10~ 3 Z & . A similar bound on 
He-shell metallicity comes from requiring a neutron den- 
sity capable of supporting a i^-induced r-process Q . As 
in this process, 9 Be synthesis requires both hard v spec- 
tra and D e o D x oscillations characteristic of an inverted 
mass hierarchy for v reactions to produce an adequate 
Y n . 

In summary, 9 Be synthesis in the He zones of mod- 
els like ull/ul5 requires an accompanying r-process and 
-E-cxpi 0.3 B. Severe fallback of inner layers was found 
to occur in all ull/ul5 models with -E CX pi ^ 0.3 B. Ta- 
ble U thus only gives ranges or upper limits for the yields 
when production partially or entirely occurs in the fall- 
back zones. The actual yields in these cases depend on 
the extent of the poorly-understood mixing and ejection 
associated with fallback. 

In the models we explored, 9 Be production is accom- 
panied by production of 7 Li and 11 B (Table [TJ) . In all 
ull/ul5 models the latter production operates as de- 
scribed in Ref. 0], through v breakup of 4 He and 12 C 
prior to shock arrival. In contrast, large fractions of 7 Li 
and n B are produced in u8.1 models by v interactions 
occurring in the O-Ne shell that has been severely al- 
tered by shock passage. Because this shell is so close 
(r ~ 1.8 x 10 8 cm) to the core, shock passage leads 
to complete disassociation of nuclei, which then are re- 
assembled into 4 He and Fe-group nuclei as the shell ex- 
pands and cools. The v irradiation of this material leads 
to 4 He( 3 He,7) 7 Be( 4 He,7) n C. The subsequent decay of 
7 Be and U C accounts for, e.g., ~70% of the 7 Li and 
^43% of the U B produced in u8.1H.l. Properties of 7 Li 
and 11 B production in all u8.1 models studied include: 1) 
weak dependence on v spectra, flavor oscillations, and the 
explosion energy and no fallback; 2) stable yields within a 
factor of ~2 of 1.6 x 10~ 7 and 8.6 x 10~ 8 M , respectively; 
3) a number-yield ratio 11 B/ 7 Li of ~0.2-0.4, comparable 
to the solar value of 0.29 18[ and distinct from the val- 
ues of ~l-3.6 found in ull/ul5 models. The last point is 
of interest because in the more massive progenitors n B 
tends to overwhelm all other yields 0, Ql| . For all the 
models in Tabled production of 6 Li and 10 B is negligible 
with 6 Li/ 7 Li < 10" 4 and 10 B/ n B - (0.5-2) x 10" 2 . 

Proposed mechanisms for Be production can be tested 
against observations. The GCR mechanism (e.g., [4|) is 
severely constrained by a recent study @ showing that 
[Be/Fe] = log (Be/Fe) -log (Be/Fe) Q - 0±0.5 for a large 
sample of stars covering [Fe/H] ~ —3.5 to —0.5. Using 
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the yields in Table Q] and the solar abundances in Ref. [l8j], 
one obtains [Be/Fc] ~ ± 0.2 for at least three cases: (I) 
all the u8.1 models with the hard v emission spectra, (II) 
ul 1*11.1 with the maximum possible Fe yield, and (III) 
ullH.l, ull*H.3, and ul5*H.3 if their actual Fe yields 
are ^10% of the upper limits. 

Our new Be production mechanisms will be most 
prominent at low metallicities, when the ISM was en- 
riched on average by a single CCSN. If the ejecta 
of a CCSN mixes with a total mass of hydrogen of 

~lO 3 (£ exp i/O.lB) 6 / 7 M [H, one finds [Fe/H] 3.5 

to -2.8 (I), -1.4 (II), and -3 to -2.4 (III) for the 
above cases. The corresponding enrichment of Li is 
A(Li) = log(Li/H) + 12 - 1-1.8. This is well below the 
level due to BBN and thus consistent with the plateau 
A(Li) ~ 2.2 observed at these metallicities [2J. Our mod- 
els also give A(B) ~ 0.6-1.3, 1.4-1.9, and 1.3-1.9 with 
B/Bc - (3-9) x 10 2 , 20-85, and 10 2 -10 3 in cases I, II, and 
III, respectively. A small number of data on B exist at 
[Fe/H] < -1.4, showing A(B) ~ 0-1.8 but typically with 
rather low B/Be values of ~10-20 [Hj]. Large B/Be val- 
ues have been observed but usually attributed to greater 
depletion of Be relative to B in stars [2l| . Simultaneous 
observations of Be and B carried out for more stars at 
[Fe/H] < —1.4 could test this interpretation, versus the 
intrinsic high yield ratio of B/Be for our models. 

In conclusion, we have identified two I'-driven CCSN- 
associated mechanisms for producing 9 Be. The first oper- 
ates in low-mass, compact progenitors, and is only mod- 
erately dependent on v spectra and flavor oscillations 



and insensitive to metallicity and the explosion energy. 
The second mechanism operates in heavier progenitors 
and requires special conditions: hard v emission spectra, 
v e <H» v x oscillations, low metallicities, and low explo- 
sion energies. The first three conditions are also those 
required for a z^-driven r-process to occur in these pro- 
genitors. Either mechanism could account for the Be seen 
in the early Galaxy. The conditions we identified for 9 Be 
synthesis are likely to arise in other astrophysical sites. 
For example, the merger of two white dwarfs could pro- 
duce a CCSN similar to the low-mass, compact progeni- 
tors studied here (e.g., (22|). In addition, significant ejec- 
tion of 9 Be was found in recent simulations of the merger 
of two neutron stars (NSs) (23[. It is quite possible that 
a variety of sources contributed to Be enrichment in the 
Galaxy, in contrast to exclusive production by GCRs. 
As Be production may be associated with both the v- 
induced r-process at low metallicities ([Fe/H] < —3) and 
the r-process in NS mergers at higher metallicities (per- 
haps [Fe/H] > —2.5), simultaneous observations of Be 
and r-process elements ought to be pursued. 
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